ESTROGEN RECEPTOR ALLELES THAT ARE PREDICTIVE OF INCREASED 
SUSCEPTIBILITY TO BONE FRACTURE 



CROSS-REFERENCES TO RELATED APPLICATIONS 
5 This application claims the benefit of U.S. Provisional Application 

No. 60/466,963, filed April 30, 2003, and U.S. Provisional Application No. 60/469,618, 
filed May 9, 2003. 

FIELD OF THE INVENTION 
The present invention relates to methods for determining susceptibility to bone 
10 fracture in a subject, methods of treating a subject to reduce the risk of bone fracture, 
methods of formulating a treatment regimen to decrease the risk of bone fracture, and kits 
for use in determining susceptibility to bone fracture. 

BACKGROUND OF THE INVENTION 
Osteoporosis is characterized by low bone mineral density (BMD), deterioration 

15 of the microarchitecture of bone and subsequent increased fracture. Twin and family 
studies have suggested that BMD has a strong genetic component, besides being 
influenced by nutritional and life-style factors. Osteoporosis is regarded as a complex 
genetic trait, which means that variants of several genes underlie the variability of the 
phenotype. Among the candidate genes in relation to BMD are the genes for collagen 

20 type IA1 (COLIA1), the Vitamin D receptor (VDR) and the estrogen receptor a (ERa). 
ERa is also known as estrogen receptor I (ESR I). Polymorphisms in the genes for the 
VDR and the ERa have been examined in relation to BMD. Although contrary reports 
have been published, two meta-analyses have shown a weak relation between the VDR 
gene and BMD. We and others have found a significant association between VDR 
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polymorphisms and fracture risk while other studies could not confirm such an 
association. Also contrary reports regarding the contribution of polymorphisms in the 
ERa gene to BMD and fracture risk have been published. A recent meta-analysis has 
shown a relation between the ERa gene and BMD and fracture risk. 
5 The vitamin D receptor gene (12ql2) comprises inherited polymorphisms 

between exon 7 and 3' UTR of the VDR gene. These alleles are denoted B/b, A/a and T/t 
for restriction enzyme sites Bsml, Apal and Taql respectively (or enzymatic or chemical 
procedures having similar specificity), where a lower case letter denotes the presence of a 
wild type restriction site which is capable of being cleaved, and a capital letter denotes 

10 the presence of a mutant restriction enzyme site which is not capable of being cleaved by 
the relevant restriction enzyme. For the purposes of the present invention, determination 
of which alleles are present in a particular gene may be referred to as determining the 
genotype of a subject for a particular gene. It is apparent from the above that each copy 
of the vitamin D receptor gene comprises a specific combination of the three alleles, this 

1 5 combination being referred to as the haplotype of the gene. For example, the haplotype 
may be baT, indicating the presence of cleavable Bsml and Apal sites, and a non- 
cleavable Taql site. Direct haplotyping of the VDR gene has allowed five different 
haplotypes to be determined, of which three are common. 

Likewise the estrogen receptor a gene (also known as estrogen receptor I) has 

20 alleles denoted P/p and X/x for restrictions sites Pvull and Xbal respectively (or 
enzymatic or chemical procedures having similar specificity), where a lower case letter 
denotes the presence of a wild type restriction site which is capable of being cleaved, and 
a capital letter denotes the presence of a mutant restriction enzyme site which is not 
capable of being cleaved by the relevant restriction enzyme. For the purpose of the 

25 present invention, determination of which alleles are present in a particular gene may be 
referred to as determining the genotype of a subject for a particular gene. It is apparent 
from the above that each copy of the estrogen receptor a gene comprises a specific 
combination of the alleles, this combination is referred to as the haplotype of the gene. 
For example, the haplotype may be px, indication the presence of cleavable Pvul and 

30 Xbal sites. Direct haplotyping of the ERa gene has allowed three different haplotypes to 
be determined. 
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VDR and ERa are interesting genes because several interactions between the 
vitamin D and estrogen endocrine system have been described. 1,25-Dihydroxyvitamin 
D 3 (l,25-(OH)2D3) and 1 7Ji-estradiol (E2) have a mutual effect on their biosynthesis and 
receptor expression. Also some genetic studies found an interaction between ERa and 
5 VDR genotypes with respect to BMD. Suarez et al. found an interactive effect of ERa 
and VDR gene polymorphisms on growth in infants. 

So far most genetic studies on osteoporosis have focused on BMD as the primary 
endpoint and not on the clinically more relevant endpoint of fractures. In contrast, the 
invention disclosed in the present application utilizes the unexpected association between 
10 specific ERa and VDR genotypes and the most typical osteoporotic fracture, the vertebral 
fracture. 

SUMMARY OF THE INVENTION 
In a first aspect, the present invention provides methods for determining 
susceptibility to bone fracture in a mammalian subject comprising at least one estrogen 

15 receptor a gene comprising a PvwII site and a Xbal site, wherein the PvwII site can exist 
as a P or p allelic form, and the Xbal site can exist as an X or x allelic form, wherein the 
method comprises analyzing nucleic acid molecules (e.g., genomic DNA) obtained from 
the mammalian subject to determine which of the P, p X, and x alleles of the estrogen 
receptor a gene are present, wherein the presence of a haplotype comprising the p and x 

20 alleles is indicative of an increased susceptibility to bone fracture. 

In some embodiments the method comprises analyzing nucleic acid molecules of 
a mammalian subject to determine the presence of the px haplotype of the estrogen 
receptor a gene, wherein the presence of said px haplotype is indicative of an increased 
susceptibility to bone fracture. In some embodiments, the method further comprises 

25 analyzing nucleic acid molecules of a mammalian subject to determine whether an allele 
of the vitamin D receptor gene is present which allele is indicative of an increased 
susceptibility to bone fracture. In some embodiments, the method comprises analyzing 
nucleic acid molecules of a mammalian subject to determine which of the B/b, A/a and 
T/t alleles of the Bsml, Apal and Taql sites of the vitamin D receptor gene are present, 

30 wherein the presence of a haplotype comprising at least of one of the baT alleles is 
indicative of an increased susceptibility to bone fracture. The method may comprise 
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determining whether more than one of the above alleles is present. The subject may be 
further classified as heterozygous or homozygous for one or more of these alleles. 
Preferably, the method comprises the additional step of determining whether the allele(s) 
present are associated with risk of bone damage, wherein presence of the b, a or T alleles 
5 is associated with increased risk of bone damage. Homozygosity for the a or b allele may 
further increase the susceptibility to bone damage in a subject, while homozygosity for 
the t allele may further decrease susceptibility. Preferably, said method comprises 
determining whether the haplotype of the subject is associated with risk of bone damage, 
wherein the haplotype baT is associated with high risk of bone damage. A subject 

10 homozygous for said haplotype may be at higher risk of bone damage than those 
heterozygous for the haplotype. 

In some embodiments, the method comprises analyzing nucleic acid molecules of 
a mammalian subject to determine the presence of the baT haplotype of the vitamin D 
receptor gene, wherein the presence of said baT haplotype is indicative of an increased 

15 susceptibility to bone fracture. In some embodiments the method further comprises 
determining the copy number of the P/p or X/x alleles of the estrogen receptor a gene 
and/or the B/b, A/a or T/t alleles of the vitamin D receptor gene. 

In other embodiments, the method comprises comparing the allele(s) present in 
the mammalian subject with genotypes of the estrogen receptor a or vitamin D receptor 

20 genes having known degrees of risk of bone fracture. 

In some embodiments the method is performed in vitro. 

In some embodiments the method is performed on a blood or tissue sample of a 

subject. 

In a second aspect of the invention there is provided methods of treating a 
25 mammalian subject to prevent or reduce the risk of bone fracture, wherein the 
mammalian subject comprises at least one estrogen receptor a gene comprising a Pvull 
site and a Xbal site, wherein the Pvull site can exist as a P or p allelic form, and the Xbal 
site can exist as an X or x allelic form, the methods each comprising analyzing nucleic 
acid molecules obtained from the mammalian subject to determine which of the P, p, X 
30 and x alleles of the Pvull and Xbal sites of the estrogen receptor a gene are present, 
wherein the presence of a haplotype comprising the p and x alleles is indicative of an 
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increased susceptibility to bone fracture, and treating the mammalian subject to reduce 
the risk of bone fracture if the subject has a haplotype comprising the p and x alleles. 

In some embodiments the method further comprises analyzing nucleic acid 
molecules of a mammalian subject to determine the presence of the px haplotype of the 
5 estrogen receptor a gene, wherein the presence of the px haplotype is indicative of an 
increased susceptibility to bone fracture. 

In some embodiments the method further comprises analyzing nucleic acid 
molecules of a mammalian subject to determine which of the B/b, A/a and T/t alleles of 
the Bsml, Apal and Taql sites of the vitamin D receptor gene are present, wherein the 
10 presence of a haplotype comprising at least one of the baT alleles is indicative of an 
increased susceptibility to bone fracture. In some embodiments the method further 
comprises analyzing nucleic acid molecules of a mammalian subject to determine the 
presence of the baT haplotype of the vitamin D receptor gene, wherein the presence of the 
baT haplotype is indicative of an increased susceptibility to bone fracture. 
15 In some embodiments the method further comprises determining the copy number 

of the P/p or X/x alleles of the estrogen receptor a gene and/or the B/b, A/a or T/t alleles 
of the vitamin D receptor gene. 

In some embodiments the method comprises comparing the alleles present in the 
mammalian subject with genotypes of the estrogen receptor a gene and/or vitamin D 
20 receptor gene having known degrees of risk of bone fracture. 

In some embodiments of the method suitable treatments include modifications to 
lifestyle, regular exercise, changes in diet or pharmaceutical preparations. 

In some embodiments of the method the mammalian subject is a human. 
In some embodiments of the method the mammalian subject is female. 
25 In a third aspect of the invention there is provided methods of formulating a 

treatment regimen to decrease the risk of bone fracture in a mammalian subject, wherein 
the methods each comprise analyzing nucleic acid molecules of a mammalian subject to 
determine whether a px haplotype of an estrogen receptor a gene is present, wherein said 
haplotype is associated with risk of bone fracture, and formulating a treatment regimen to 
30 decrease the risk of bone fracture in the mammalian subject. 

In some embodiments the method further comprises analyzing nucleic acid 
molecules of a mammalian subject to determine the presence of the baT haplotype of the 
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vitamin D receptor gene, wherein said haplotype is associated with risk of bone fracture, 
and formulating a treatment regimen to decrease the risk of bone fracture based on said 
haplotype. 

In some embodiments the method further comprises comparing the alleles present 
5 in the mammalian subject with genotypes of the estrogen receptor a gene and/or 
vitamin D receptor gene having known degrees of risk of bone fracture. 

In some embodiments the method comprises administering an appropriate 
treatment, such as modifications to lifestyle, regular exercise, changes in diet or 
pharmaceutical preparations. 

10 In a fourth aspect of the invention there is provided methods of determining 

susceptibility to bone fracture in a mammalian subject, wherein the methods each 
comprise the step of utilizing a kit to determine whether a px haplotype of an estrogen 
receptor a gene is present in a mammalian subject, wherein said kit comprises (i) one or 
more nucleic acid primer molecules for amplification of a portion of the estrogen receptor 

15 a gene, (ii) means for determining whether the px haplotype of said gene is present, and 
(iii) means for indicating a correlation between the presence of the px haplotype and risk 
of bone fracture, and wherein the presence of the px haplotype in the mammalian subject 
is indicative of susceptibility to bone fracture. 

The association between the presence of the px haplotype and susceptibility to 

20 bone fracture is independent of BMD. In some embodiments the method further 
comprises the step of determining whether the baT haplotype of the vitamin D receptor 
gene is present in the mammalian subject, said kit further comprising (i) one or more 
nucleic acid primer molecules for amplification of a portion of the vitamin D receptor 
gene, (ii) means for determining whether the baT haplotype of the vitamin D receptor 

25 gene is present and (iii) means for indicating a correlation between the presence of the 
baT haplotype and risk of bone fracture. 

In a fifth aspect of the invention there is provided a kit for determining 
susceptibility to bone fracture in a subject, said kit comprising (i) one or more nucleic 
acid primer molecules for amplification of a portion of an estrogen receptor a gene, 

30 (ii) means for determining whether a px haplotype of said gene is present; and (iii) means 
for indicating a correlation between the presence of said haplotype and risk of bone 
fracture. 
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In some embodiments the kit further comprises (i)one or more nucleic acid 
primer molecules for amplification of a portion of a vitamin D receptor gene, (ii) means 
for determining whether a baT haplotype of the vitamin D receptor gene is present, and 
(iii) means for indicating a correlation between the presence of the baT haplotype and risk 
5 of bone fracture. 

In some embodiments of the first to fifth aspects the haplotype of the estrogen 
receptor a gene is determined by amplification of a portion of the first intron of the 
estrogen receptor a gene, followed by restriction enzyme digestion; or any other 
technique suitable for determining the genotype of the subject. 
10 In some embodiments of the first to fifth aspects the haplotype of the vitamin D 

receptor is determined by amplification of a portion of the vitamin D receptor gene 
between exon 7 and the 3' untranslated region (UTR), followed by restriction enzyme 
digestion; or any other technique suitable for determining the genotype of the subject. 

BRIEF DESCRIPTION OF THE DRAWINGS 
15 The foregoing aspects and many of the attendant advantages of this invention will 

become more readily appreciated as the same become better understood by reference to 
the following detailed description, when taken in conjunction with the accompanying 
tables and drawings, wherein: 

FIGURE 1 shows VDR gene: direct molecular haplotyping. For details see 
20 Example 1. 

FIGURE 2 shows ER gene: direct molecular haplotyping. For details see 
Example 1. 

FIGURE 3 shows odds ratios and numbers of vertebral fractures according to ERa 
haplotype 1 (px) and VDR haplotype 1 (baT) genotype. 'Reference' includes ERa and 
25 VDR genotypes 22, 23, and 33; 'hetero* includes 12 and 13 genotypes; 'homo 1 includes 
1 1 genotype. **/> value < 0.001 for allele dose association of ERa haplotype 1. */?-value 
= 0.06 for allele dose association of VDR haplotype 1. 

FIGURE 4 shows odds ratios compared to the 'double reference' group and 
numbers of vertebral fractures according to combined ERa and VDR genotypes. 
30 'Reference' includes ERa and VDR genotypes 22, 23, and 33; 'heterozygotes' include 12 
and 13; 'homozygote' includes 11. value < 0.001 for allele dose association of ERa 
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haplotype 1 (px) in VDR haplotype 1 (baT) homozygous carriers, p = 0.01 for the 
interaction term. 

FIGURE 5 shows the nucleic acid sequence of part of intron I/exdn 2 of the 
human Estrogen receptor a gene containing Pvull and xbal polymorphisms. Nucleic acid 
5 sequences to which PCR primers anneal are underlined. The nucleic acid sequence of 
this portion of the estrogen receptor a gene is also shown in SEQ ID NO: 1 . 

FIGURE 6 shows the nucleic acid sequence of the 3 ' part of the human VDR gene 
containing Bsml, Apal, and Taq\ RFLP's polymorphic sites. The nucleic acid sequences 
to which PCR primers anneal are underlined. The nucleic acid sequence of this portion of 
10 the human VDR gene is also shown in SEQ ID NO:2. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
Unless specifically defined herein, all terms used herein have the same meaning 
as they would to one skilled in the art of the present invention. Practitioners are 
particularly directed to Sambrook et al. (1989) Molecular Cloning: A Laboratory 

15 Manual, 2 nd ed., Cold Spring Harbor Press, Plainsview, New York (1989), and 
Ausubel et al., Current Protocols in Molecular Biology (Supplement 47), John Wiley & 
Sons, New York (1999), for definitions and terms of the art. 

As used herein, the following terms have the meanings defined below: Estrogen 
receptor a (abbreviated as ERa) is a nuclear hormone receptor protein that binds estrogen 

20 and is thereby activated and binds to specific sequences of nuclear DNA which serve as 
on-off switches for transcription within the cell nucleus. Representative ERa genes that 
are analyzed in the practice of the present invention hybridize under conditions of 1.0 X 
SSC, 0.5% SDS, at 30°C for 30 minutes to the portion of the representative ERa gene set 
forth in SEQ ID NO:l. Some ERa genes that are analyzed in the practice of the present 

25 invention hybridize under conditions of 1.0 X SSC, 0.5% SDS, at 45°C for 30 minutes to 
the portion of the representative ERa gene set forth in SEQ ID NO:l. Some ERa genes 
that are analyzed in the practice of the present invention hybridize under conditions of 
1.0 X SSC, 0.5% SDS, at 60°C for 30 minutes to the portion of the representative ERa 
gene set forth in SEQ ID NO:l. The ability to hybridize under the foregoing conditions 

30 can be determined, for example, by first hybridizing the nucleic acid molecules under low 
stringency conditions (e.g., 6 X SSC, 0.5% SDS, at 25°C for 12 hours), then incubating 
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the hybridized nucleic acid molecules under the stated conditions (e.g., 1.0 X SSC, 
0.5% SDS, at 30°C for 30 minutes). 

There is a polymorphism at position 397 of the nucleotide sequence (numbered 
from the transcriptional start site) of the ERa gene. The polymorphism consists of a 
5 substitution of a C residue at position 397 of the sequence with a T residue. The genetic 
polymorphism is denoted "C397T" where the number refers to the position of the 
polymorphism with respect to the nucleotide sequence; the "C" is the nucleotide residue 
present in the reference or wild type sequence; and the "T" is the nucleotide residue 
present at that position in the variant sequence. There is a polymorphism at position 351 

10 of the nucleotide sequence of the ERa gene. The polymorphism consists of a substitution 
of a G residue at position 351 for the sequence provided in SEQ ID NO:l with a A 
residue. The genetic polymorphism is denoted "G351A" where the number refers to the 
position of the polymorphism with respect to the nucleotide sequence; the M G" is the 
nucleotide residue present in the reference or wild type sequence; and the "A" is the 

15 nucleotide residue present at that position in the variant sequence. 

As used herein, the term "VDR gene" refers to a gene coding for the Vitamin D 
receptor. 

As used herein, the term "bone damage" refers to any form of structural damage 
including fractures, break or chips. The term may also include biological degradation or 

20 deterioration of bone. Typically, the term "bone damage" does not include low bone 
mineral density. This is in line with the finding that risk of bone damage is independent 
of bone mineral density. Fracture may be defined as the clinically most important 
endpoint, and thus the method of the first aspect of the invention preferably relates to a 
method for determining risk of fracture independent of bone mineral density. Although 

25 such bone damage will usually be the result of osteoporosis, it is irrelevant for the 
purposes of the present invention whether a subject has first been diagnosed as having 
osteoporosis or low bone mineral density. 

As used herein, the terms "risk of bone damage" and "susceptibility to bone 
damage" are used interchangeably and such risk or susceptibly is independent of bone 

30 mineral density. The present invention enables the identification of those individuals 
susceptible to bone fracture, and the development of therapeutic or preventative 
measures. For example, those at risk may avoid damage by modifying their lifestyle and 
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implementing bone strengthening measures, such as regular exercise and a healthy diet, 
or by taking medicaments which reduce the risk of damage. Exemplary, art-recognized 
strategies for the prevention and/or treatment of individuals who are at increased risk of 
bone fracture are described in: DeSantis, A. and A. Buchman, "Current and Emerging 
5 Therapies in Osteoporosis," Expert Opin. Pharmacother. 5(7):835-843, 2002; 
Whitfield, J.A., et al., "Bone Growth Stimulators. New Tools for Treating Bone Loss and 
Mending Fractures," Vitam. Horm. (55:1-80, 2002; Lane, J.M., et al., "Medical 
Management of Osteoporosis," Institute for Biological Sciences, National Research 
Council of Canada, Ottawa, Ontario Canada, Instr. Course Lect 52:785-789, 2003; 

10 Kim, L, "Osteoporosis in Aging-Challenges and Opportunities," J. Okla. State Med. 
Assoc. 96(3): 143-146, 2003; Body, J.J., "Management of Primary Osteoporosis," Acta 
Clin. Belg. 57(5):277-283, 2002. 

In one aspect, the present invention provides methods for determining 
susceptibility to bone fracture in a mammalian subject comprising at least one estrogen 

15 receptor a gene comprising a Pvull site and a Xbal site, wherein the Pvull site can exist 
as a P or p allelic form, and the Xbal site can exist as an X or x allelic form, wherein the 
method comprises analyzing nucleic acid molecules (e.g., genomic DNA) obtained from 
the mammalian subject to determine which of the P, p X, and x alleles of the estrogen 
receptor a gene are present, wherein the presence of a haplotype comprising the p and x 

20 alleles is indicative of an increased susceptibility to bone fracture. 

In some embodiments, the method further comprises analyzing nucleic acid 
molecules of a mammalian subject to determine whether an allele of the vitamin D 
receptor gene is present which allele is indicative of an increased susceptibility to bone 
fracture. In some embodiments, the method comprises analyzing nucleic acid molecules 

25 of a mammalian subject to determine which of the B/b, A/a and T/t alleles of the BsmI, 
Apal and Taql sites of the vitamin D receptor gene are present, wherein the presence of a 
haplotype comprising at least of one of the baT alleles is indicative of an increased 
susceptibility to bone fracture. 

The determination that an allele of a polymorphism in an estrogen receptor a gene 

30 and/or a Vitamin D receptor gene is associated with an increased or decreased 
susceptibility to bone damage may be performed using standard statistical analyses in a 
population of subjects, as described, for example in the "statistical analysis" section. An 
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exemplary method for determining the presence of one or more alleles associated with an 
increased or a decreased susceptibility to bone damage is present in a subject is described 
in detail below. An exemplary haplotype is the ERa haplotype 1 (px) of the estrogen 
receptor a and the VDR haplotype 1 (baT) of the VDR receptor. The estrogen receptor a 
5 gene comprises an inherited polymorphism at positions 397 and 351 of a nucleotide 
sequence of the ERa. The sequence of the estrogen receptor a gene and methodology on 
how to identify the C397T and G351A polymorphisms have been previously described. 

The present invention is based upon the surprising observation of a correlation 
between the presence of the T allele at position 397 and/or the A allele at position 351 of 

10 the estrogen receptor a gene and the susceptibility to/or risk of bone fracture in those 
subjects having the polymorphisms independent of BMD. A subject having the estrogen 
receptor a T allele and the A allele in addition to the VDR haplotype 1 (baT) shows a 
higher risk of fracture compared to a subject having the estrogen receptor a C and 
G alleles and a haplotype of 2 (BAt) or 3 (bAT) for the VDR polymorphism. 

15 These results could not be predicted as previous studies have shown that the 

fracture risk is dependent on bone mineral density. This fact is borne out by the results 
presented herein, where those individuals having the highest risk of bone fracture are 
those not having low bone mineral density. By screening for the alleles of the estrogen 
receptor a gene, and/or the VDR gene susceptibility to bone damage may be assessed 

20 without the need for an analysis of bone mineral density. 

In some embodiments there are provided methods for determining susceptibility 
to bone damage in a subject comprising determining the presence of T allele of the 
C397T polymorphism of the estrogen receptor a gene and/or the A allele of the G351A 
polymorphism of the estrogen receptor a gene. This fact is confirmed by the result 

25 presented herein, that those individuals at highest risk of bone damage have both 
polymorphisms. Thus, some embodiments provide methods for determining 
susceptibility to bone damage in a subject, comprising determining the presence of the 
T allele of the C397T polymorphism of the estrogen receptor a gene and the A allele of 
the G351A polymorphism of the estrogen receptor a gene, wherein the presence of at 

30 least one copy of the T allele of the C397T polymorphism and at least one copy of the 
A allele of the G351A polymorphism of the estrogen receptor a gene in the subject is 
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indicative of an increased susceptibility to bone fracture independent of bone mineral 
density. 

Typically, the method of the first aspect of the present invention comprises 
analysis of polymorphisms in estrogen receptor a gene to determine susceptibility to bone 
5 damage. The method may include analysis of DNA, and/or RNA, to determine which 
allele of a polymorphism is present. The method may include determining whether one 
or more particular alleles are present. The method may further comprise determining 
whether subjects are homozygous or heterozygous for alleles of the estrogen receptor a 
gene. 

10 Preferably, the method of the first aspect of the present invention further 

comprises analyzing nucleic acid molecules of a mammalian subject to determine 
whether an allele of the Vitamin D receptor is present which allele is indicative of an 
increased susceptibility to bone fracture. This may be performed by comparing the 
alleles present in a subject with those known to be associated with risk of bone damage. 

15 For example, a visual aid detailing alleles and the relative risk of bone damage associated 
therewith may be used to determine whether the genotype or phenotype of the subject is 
associated with a high or low risk of bone damage. 

The methods of the present invention may be performed in vitro. Preferably, the 
method is performed on a tissue or fluid sample removed from the body of the subject. 

20 Thus, the present invention relates to a non-invasive diagnostic method, the results of 
which provide an indication of susceptibility to bone fracture but do not lead to a 
diagnosis upon which an immediate medical decision regarding treatment can be made. 

. The present invention may be performed on any living subject for which it is 
desirable to determine risk of bone fracture. Most preferably, the subject is a human, 

25 preferably a female. 

The present invention may be performed using any suitable method known in the 
art. Preferably, a tissue or fluid sample is first removed from a subject. Examples of 
suitable samples include blood, mouth or cheek cells, and hair samples containing roots. 
Other suitable samples would be known to the person skilled in the art. The nucleic acid 

30 molecules are extracted from the sample, using any suitable method. For example the 
nucleic acid molecules may be extracted using the techniques described in Sambrook et al 
(Molecular Cloning - A Laboratory Manual, Cold Spring Harbor Laboratory Press). 
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Determination of the genotype or phenotype of a subject may then be carried out using 
the extracted nucleic acid molecules, employing any suitable technique, including, for 
example, Southern blot analysis followed by restriction enzyme digestion; PCR 
amplification followed by restriction enzyme digestion and, optionally, separation of 
5 digestion products by gel electrophoresis; sequencing of a relevant gene fragment by any 
suitable method; visualization of heteroduplex patterns, for example on PAA or agarose 
gels, where different patterns may indicate the presence of one or more specific alleles; 
separation of DNA fragments using denaturing gradient gels, wherein the degree of 
separation will depend upon the presence or absence of one or more polymorphic 

10 restriction sites; separation using SSCP analysis, the patterns of which will depend upon 
the presence or absence of one or more polymorphic restriction sites, use of allele specific 
oligonucleotides, hybridization patterns of which some will be specific for various 
combinations of alleles; methods such as OLA, Taqman or dot-bloting for the detection 
of known mutations; visualization of DNA sites using fluorescent labeled probes for 

15 alleles of interest; and RFLP analysis. Other suitable methods will be known to persons 
skilled in the art. 

Where it is desirable to use particular restriction enzymes in performing the 
present invention, the skilled person will understand that enzymatic or chemical 
procedures having similar specificities may also be used. For example, restriction 
20 enzymes having similar specificity (isoschizomers) to those described herein may be 
used, or chemical degradation procedures with DNA or RNA cutting specificity. 

Other techniques suitable for determining the genotype or phenotype of a subject 
may be used in the present invention. 

Amplification is preferably carried out by polymerase chain reaction (PCR) 
25 techniques, to produce copies of the amplified region, where the fragment is of the 
estrogen receptor a gene, the fragments are at least about 20, preferably at least about 50, 
about 70, about 100, about 150, or about 200 bases in length. Where the fragment to be 
amplified is of the VDR gene, PCR primers may be selected to amplify a fragment which 
is at least about 20, preferably at least about 50, about 70, about 100, about 150, or about 
30 200 bases in length. 

Exemplary PCR primers are at least about 10 nucleotides in length, preferably at 
least about 1 5 nucleotides or at least 20 nucleotides in length, and are complementary to 
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any stretch of at least about 10 nucleotides of the sequence to be amplified. PCR 
techniques are well know in the art, and it is within the ambit of the skilled person to 
identify primers for amplification of the appropriate region of the above genes, namely 
the first intron of the estrogen receptor a gene and the region from exon 7 to the 3 f UTR 
5 of the VDR gene. A preferred technique is single base extension, and for this method it is 
only necessary to amplify a fragment including the polymorphic site. Thus, amplification 
of the region immediately surrounding the C397T and G351A polymorphic site is 
necessary. Exemplary PCR techniques are described in EP-A-0200362 and EP-A- 
0201184. 

10 In a preferred feature of the first aspect, there is provided a method for 

determining susceptibility to bone fracture in a subject, said method comprising 
amplifying a fragment comprising a portion of the first intron of the estrogen receptor a 
gene, and determining which allele(s) in the estrogen receptor a gene is/are present. 
Primers suitable for amplification of said portion of the estrogen receptor a gene would 
15 be readily available to a person skilled in the art. Examples for such primers include: 
SEQ ID NO:3 : 5'-GATATCCAGGGTTATGTGGCA-3' 
SEQ ID NO:4: 5 1 - AGGTGTTGCCTATTAT ATT AACCTTG A-3 ' 

To determine which allele of the VDR gene is present, at least a portion of the 
region from exon 7 to 3' UTR of the VDR may be amplified followed by the 
20 determination of the presence of the Bsml, Apal and Taql sites. Suitable primers include: 
SEQ ID NO:5 5 '-C AACC AAGACTAC AAGTACCGCGTC AGTGA-3 ' 

and/or 

SEQ ID NO:6 5 '-TTTGGCTCC AATC AGAT AC ATGGGA-3 ' 

SEQ ID NO:7 5 f -GCAACTCCTCATGGCTGAGGTCTC-3 f 

25 In a second aspect of the invention there is provided methods of treating a 

mammalian subject to prevent or reduce the risk of bone fracture, wherein the 
mammalian subject comprises at least one estrogen receptor a gene comprising a Pvull 
site and a Xbal site, wherein the Pvull site can exist as a P or p allelic form, and the Xbal 
site can exist as an X or x allelic form, the methods each comprising analyzing nucleic 

30 acid molecules obtained from the mammalian subject to determine which of the P, p, X 
and x alleles of the Pvull and Xbal sites of the estrogen receptor a gene are present, 
wherein the presence of a haplotype comprising the p and x alleles is indicative of an 
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increased susceptibility to bone fracture, and treating the mammalian subject to reduce 
the risk of bone fracture if the subject has a haplotype comprising the p and x alleles. 

In some embodiments the method further comprises analyzing nucleic acid 
molecules of the mammalian subject to determine which of the B/b, A/a and T/t alleles of 
5 the BsmI, Apal and Taql sites of the vitamin D receptor gene are present, wherein the 
presence of a haplotype comprising at least one of the baT alleles is indicative of an 
increased susceptibility to bone fracture. In some embodiments the method further 
comprises analyzing nucleic acid molecules of the mammalian subject to determine the 
presence of the baT haplotype of the vitamin D receptor gene, wherein the presence of the 
10 baT haplotype is indicative of an increased susceptibility to bone fracture. 

In some embodiments the subject is diagnosed as being at risk of bone fracture, 
preferably using the method of the first aspect of the invention. In this aspect the 
reduction in risk of bone fracture includes any means of reducing risk of bone fracture in 
a subject. 

15 Therapy may be in the form of preventative or palliative care. A preferred 

method of treatment is prescribing or administering an agent that reduces the 
susceptibility of a subject to bone fracture. For example, the treatment may comprise 
prescribing or administering Fosamax® (sodium alendronate) or Fosteo® (parathyroid 
hormone) or anabolic steroids, bisphosphonates, vitamin D preparations or other 

20 treatments known to those skilled in the art. Other suitable treatment which may be 
prescribed or administered alongside Fosamax or Fosteo, or as an alternative thereto, 
treatment also includes modifications to lifestyle, regular exercise and changes in diet to 
strengthen bones, and hormone therapy. Other suitable treatments including 
pharmaceutical preparations to reduce bone loss as known to physicians and persons 

25 skilled in the art. Examples include anabolic steroids, bisphosphonates, vitamin D 
preparations, calcium supplements and Hormone Replacement Therapy. 

Administration of the medicament is accomplished by any effective route, 
e.g., orally or parenteral. Methods for parenteral delivery include topical, intra-arterial, 
subcutaneous, intramedullary, intravenous, or intranasal administration. Oral 

30 administration followed by subcutaneous injection would be the preferred routes of 
uptake, also long acting immobilizations would be used. In addition to the active 
ingredients, these medicaments may contain suitable pharmaceutically acceptable carriers 
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comprising excipients and other compounds that facilitate processing of the active 
compounds into preparations which can be used pharmaceutically. Further details of 
techniques for formulation and administration may be found in the latest edition of 
"Remington's Pharmaceutical Sciences" (Mack Publishing Co, Easton PA). 
5 Medicaments for oral administration can be formulated using pharmaceutically 

acceptable carriers well known in the art, in dosages suitable for oral administration. 
Such carriers enable the pharmaceutical compositions to be formulated as tablets, pills, 
dragees, capsules, liquids, gels, syrups, slurries, suspensions, etc., suitable for ingestion 
by the patient. 

10 Medicaments for oral use can be obtained through combination of active 

compounds with solid excipient, optionally grinding a resulting mixture, and processing 
the mixture of granules, after adding suitable addition compounds, if desired, to obtain 
tablets or dragee cores. Suitable excipients are carbohydrate or protein fillers. These 
include, but are not limited to sugars, including lactose, sucrose, mannitol, or sorbitol, 

1 5 starch from corn, wheat, rice, potato, or other plants; cellulose such as methyl cellulose, 
hydroxypropylmethyl-cellulose, or sodium carboxymethylcellulose; and gums including 
arabic and tragacanth; as well as proteins, such as gelatin or collagen. If desired, 
disintegrating or solubilizing agents may be added, such as the cross-linked polyvinyl 
pyrrolidone, agar, alginic acid, or a salt thereof, such as sodium alginate. 

20 Dragee cores are provided with suitable coating such as concentrated sugar 

solutions, which may also contain gum arabic, talc, polyvinylpyrrolidone, carbopol gel, 
polyethylene, glycol, and/or titanium dioxide, lacquer solutions, and suitable organic 
solvents or solvent mixtures. Dyestuffs or pigments may be added to the tablets or dragee 
coatings for product identification or to characterize the quantity of active compound 

25 (i.e., dosage). 

Medicaments, which can be used orally, include push-fit capsules made of gelatin, 
as well as soft, sealed capsules made of gelatin and a coating such as glycerol or sorbitol. 
Push-fit capsules can contain active ingredients mixed with filler or binders such as 
lactose or starches, lubricants such as talc or magnesium stearate, and, optionally, 
30 stabilizers. In soft capsules, the active compounds may be dissolved or suspended in 
suitable liquids, such as fatty oils, liquid paraffin, or liquid polyethylene glycol with or 
without stabilizers. 
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Medicaments for parenteral administration include aqueous solutions of active 
compounds. For injection, the medicaments of the invention may be formulated in 
aqueous solutions, preferably in physiologically compatible buffers such as Hank's 
solution, Ringer's solution, or physiologically buffered saline. Aqueous injection 
5 suspensions may contain substances, which increase the viscosity of the suspension, such 
as sodium carboxymethyl cellulose, sorbitol, or dextran. Additionally, suspensions of the 
active compounds may be prepared as appropriate oily injection suspensions. Suitable 
lipophilic solvents or vehicles include fatty oils such as sesame oil, or synthetic fatty acid 
esters, such as ethyl oleate or triglycerides, or lipsomes. Optionally, the suspension may 
10 also contain suitable stabilizers or agents, which increase the solubility of the compounds 
to allow for the preparation of highly concentrated solutions. 

For topical or nasal administration, penetrants appropriate to the particular barrier 
to be permeated are used in the formulation. Such penetrants are generally known in the 
art. 

15 Medicaments useful for treating mammalian subjects susceptible to bone fracture 

may be manufactured in a manner similar to that know in the art (e.g., by means of 
conventional mixing, dissolving, granulation, dragee-making, levigating, emulsifying, 
encapsulation, entrapping or lyophilizing processes). The medicaments may also be 
modified to provide appropriate releases characteristics, e.g., sustained release or targeted 

20 release, by conventional means e.g., coating. 

The medicaments may be provided as a salt and can be formed with many acids, 
including but not limited to hydrochloric, sulfuric, acetic, lactic, tartaric, malic, succinic, 
etc. Salts tend to be more soluble in aqueous or other protonic solvents that are the 
corresponding free base forms. In other cases, the preferred preparation may be a 

25 lyophilized powder in 1 mM-50mM histidine, 0.1%-2% sucrose, 2%-7% mannitol at a 
pH range of 4.5 to 5.5, that is combined with buffer prior to use. 

After such medicaments formulated in an acceptable carrier have been prepared, 
they can be placed in an appropriate container and labeled for treatment of an indicated 
condition. 

30 Medicaments suitable for prevention or reduction of bone damage include 

compositions wherein the active ingredients are contained in an effective amount to 
achieve the intended purpose. The amount actually administered will be dependent upon 
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the individual to which treatment is to be applied, and will preferably be an optimized 
amount such that the desired effect is achieved without significant side-effects. The 
determination of a therapeutically effective dose is well within the capability of those 
skilled in the art. Of course, the skilled person will realize that divided and partial doses 
5 are also within the scope of the invention. 

For any compound, the therapeutically effective dose can be estimated initially 
either in cell culture assays or in any appropriate animal model (e.g., primate, rats and 
guinea pigs for hypertension and other small laboratory animals). These assays should 
take into account receptor activity as well as downstream processing activity. The animal 

10 model is also used to achieve a desirable concentration range and route of administration. 
Such information can then be used to determine useful doses and routes for 
administration in humans. 

A therapeutically effective amount refers to that amount of agent, which 
ameliorates the symptoms or condition. Therapeutic efficacy and toxicity of such 

15 compounds can be determined by standard pharmaceutical procedures, in cell cultures or 
experimental animals (e.g., EDso, the dose therapeutically effective in 50% of the 
population; and LD 5 o> the dose lethal to 50% of the population). The dose ratio between 
therapeutic and toxic effects is the therapeutic index, and it can be expressed as the ratio 
ED 5 o/LD 50 . Medicaments, which exhibit large therapeutic indices, are preferred. The 

20 data obtained from cell culture assays and animal studies is used in formulating a range of 
dosage for human use. The dosage of such compounds lies preferably within a range of 
circulation concentrations that include the ED50 with little or no toxicity. The dosage 
varies within this range depending upon the dosage form employed, sensitivity of the 
patient, and the route of administration. 

25 The exact dosage is chosen by the individual physician in view of the patient to be 

treated. Dosage and administration are adjusted to provide sufficient levels of the active 
moiety or to maintain the desired effect. Long acting medicaments might be administered 
every 3 to 4 days, every week, or once every two weeks depending on half-life and 
clearance rate of the particular formulation. Guidance as to particular dosages and 

30 methods for delivery is provided in the literature (see U.S. Patent Nos. 4,657,760; 
5,206,344 and 5,225,212, herein incorporated by reference). 
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In a third aspect of the present invention, there is provided a method of 
formulating a treatment regimen to decrease the risk of bone fracture. In some 
embodiments the method comprises analyzing nucleic acid molecules of a mammalian 
subject to determine the presence of the px haplotype of the ERa (haplotype 1) and the 
5 baT haplotype of the VDR gene (haplotype 1). This may be done according to the 
methods provided in the first aspect of the invention. 

In a third aspect of the present invention, there is provided a method of 
formulating a treatment regimen to decrease the risk of bone fracture. In some 
embodiments the method comprises analyzing nucleic acid molecules and/or proteins of a 

10 mammalian subject to determine the presence of the px haplotype of the estrogen receptor 
a gene, wherein said haplotype is associated with risk of bone fracture, and formulating a 
treatment regimen to decrease the risk of bone fracture. The nucleic acid molecules can 
be analyzed using, for example, the techniques described in connection with the methods 
provided in the first aspect of the invention. 

15 The effect of a therapeutic or preventative agent may depend on the underlying 

cause of the bone damage, in some cases it may be preferable to avoid the use of certain 
treatments for example, the presence or absence of particular alleles of a gene will 
provide a useful indication as to which is the most appropriate preventative measure. 
This aspect of the present invention may also be useful for identifying agents which may 

20 be used in the treatment of bone damage. 

Thus, in some embodiments, the present invention provides methods for 
predicting the response of a mammalian subject to treatment to decrease the risk of bone 
damage that comprises the steps of: (1) determining which allele(s) of polymorphisms in 
an ERa gene, or the combination of an ERa gene and a VDR gene a mammalian subject 

25 has; and (2) predicting the response of the subject to treatment, wherein the absence of an 
allele of a polymorphism in an ERa gene or a VDR gene that is associated with an 
increased risk of bone damage is indicative that treatment of the mammalian subject is 
unlikely to reduce the risk of bone damage. The polymorphisms in the ERa gene that is 
associated with an increased risk of bone damage may be the C397T and or G351A 

30 polymorphisms. The polymorphisms in the VDR gene that are associated with an 
increased risk of bone damage may be the baT polymorphisms. 
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In some embodiments, the method comprises the steps of: (1) determining that a 
mammalian subject has an increased susceptibility to bone damage by using a method 
comprising the step of determining the presence in the subject of one or more allele(s) of 
polymorphisms in an ERa gene and a VDR gene that are associated with increased 
5 susceptibility to bone damage; and (2) predicting the response of the subject to treatment 
wherein the absence of an allele of a polymorphism in ERa that is associated with an 
increased risk of bone damage is indicative that treatment of the subject is unlikely to 
reduce the risk of bone damage. The polymorphism in the ERa gene that is associated 
with an increased risk of bone damage may be the C397T and/or the G351A 

10 polymorphisms. The polymorphism in the VDR gene that is associated with an increased 
risk of bone damage may be the baT polymorphism. 

In a fourth aspect of the invention there is provided methods of determining 
susceptibility to bone fracture in a mammalian subject, wherein each method comprises 
the step of utilizing a kit to determine whether a px haplotype of an estrogen receptor a 

15 gene is present in a mammalian subject, wherein said kit comprises (i) one or more 
nucleic acid primer molecules for amplification of a portion of the estrogen receptor a 
gene, (ii) means for determining whether the px haplotype of said gene is present, and 
(iii) means for indicating a correlation between the presence of the px haplotype and risk 
of bone fracture, and wherein the presence of the px haplotype in the mammalian subject 

20 is indicative of susceptibility to bone fracture. 

Some embodiments of the methods further comprise the step of determining 
whether a baT haplotype of a vitamin D receptor gene is present in the mammalian 
subject, said kit further comprising (i) one or more nucleic acid primer molecules for 
amplification of a portion of the vitamin D receptor gene, (ii) means for determining 

25 whether the baT haplotype of the vitamin D receptor gene is present, and (iii) means for 
indicating a correlation between the presence of the baT haplotype and risk of bone 
fracture. 

Details of compounds of the Kit include those described below in relation to the 
fifth aspect. 

30 The Method of the fourth aspect may be performed using any suitable methods 

known in the art, for example those as discussed above in relation to the method of the 
first aspect of the invention. 
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In a fifth aspect of the present invention, there is provided kits for use in 
determining which allele(s) of a polymorphism of an ERa gene or the combination of an 
ERa gene and a VDR gene are present in a mammalian subject. The kits include (i) one 
or more nucleic acid primer molecules for amplification of a portion of an estrogen 
5 receptor a gene, (ii) means for determining whether a px haplotype of said gene is 
present; and (iii) means for indicating a correlation between the presence of said 
haplotype and risk of bone fracture. In some embodiments the kits include (i) one or 
more nucleic acid primer molecules for amplification of a portion of a vitamin D receptor 
gene, (ii) means for determining whether a baT haplotype of the vitamin D receptor gene 

10 is present, and (iii) means for indicating a correlation between the presence of the baT 
haplotype and risk of bone fracture. 

In some embodiments, the kit provides a method for determining susceptibility to 
bone damage by determining the presence of alleles of ERa or alleles of ERa and VDR 
that are associated with an increased or decreased susceptibility to bone damage. In 

15 further embodiments, the kit may be used in a method of predicting the response of a 
subject to treatment by determining the presence of alleles of ERa or alleles of ERa and 
VDR that are associated with an increased or decreased susceptibility to bone damage. 

Preferably, the kit also comprises means for indicating a correlation between the 
allele(s) and risk of bone damage, or for predicting the response of a mammalian subject 

20 to treatment. Thus, the kit may contain written indicia providing information to the user 
for interpreting the results of the analysis. For example the written indicia may explain 
that the presence of the Pvull and Xbal restriction sites at the p and x polymorphisms of 
the ERa gene in a subject indicates the presence of the T and A alleles at these 
polymorphic sites. 

25 Preferably, the primer molecules are suitable for amplification of at least a portion 

of the ERa gene or the ERa and VDR genes. Primers suitable for amplification of a 
portion of the ERa gene would be readily available to a person skilled in the art. 
Examples of suitable primers are described above. For example, suitable primers may 
amplify a fragment comprising a portion of the ERa gene, and determine which allele(s) 

30 of the C397T and G351A (px) polymorphisms in ERa is/are present. Examples of such 
primers include: 

SEQ ID NO:3 : S'-GATATCCAGGGTTATGTGGCA-S' 
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SEQ ID NO:4: 5' -AGGTGTTGCCTATTATATTAACCTTGA-3 1 

Primers suitable for amplification of a portion of the VDR gene would also be 
readily available to a person skilled in the art. For example, suitable primers may amplify 
a portion of the VDR gene to determine which allele(s) of the baT polymorphisms of the 
5 VDR gene is present. Suitable primers include: 

SEQ ID NO:5 5'-CAACCAAGACTACAAGTACCGCGTCAGTGA-3' 

and/or 

SEQ ID NO:6 5 '-TTTGGCTCC AATC AGATAC ATGGGA-3 1 

SEQ ID NO:7 5'-GCAACTCCTCATGGCTGAGGTCTC-3 , 

10 Means for determining which allele(s) is/are present in the ERa gene, and/or VDR 

gene may include any reagents or molecules necessary for use in any of the methods 
described above. For example, where PCR followed by DNA digestion is used, said 
means preferably include PCR reagents and one or more of the Pvull and/or Xbal 
restriction enzymes. Where the method employs Southern Blotting, heteroduplex 

15 visualization, or fluorescent labeling techniques for example, probes which bind to the 
appropriate region of the ERa gene, and/or VDR gene may be included. Where 
necessary, such probes may be labeled to allow detection, for example by nick- 
translation, radio- or fluorescent-labeling, or random primer extension whereby the non- 
labeled nucleotides serve as a template for the synthesis of labeled molecules. Other 

20 methods for labeling probes are well known in the art. 

The means for correlating the allele present with risk of bone damage, or for 
predicting response to treatment, may be in the form of a chart or visual aid. The chart or 
visual aid may indicate that presence of the p and x alleles of the ERa gene and the baT 
alleles of the VDR gene is associated with increased risk of bone damage. In a preferred 

25 feature of the fifth aspect, the kit may also comprise control DNA samples, for 
comparison with DNA sequences of a subject. The control samples may comprise the 
sequence of one or more alleles of the ERa and/or VDR genes. 

The following examples are provided for the purpose of illustration, not limiting, 
the present invention. 
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EXAMPLE 1 

This Example describes the identification of polymorphisms in the ERa gene and 
the VDR gene that are positively correlated with increased susceptibility to bone fracture 
in human beings. 

5 STUDY SUBJECTS 

All women included in this study were part of a population-based cohort study of 
subjects aged 55 years or more, living in the Ommoord district of the city of Rotterdam in 
The Netherlands. The objective of the study was to document the occurrence of disease 
in the elderly in relation to several potential determinants (1). A total of 10,275 persons 

10 were invited for baseline examination in 1990. Of those 7,983 (61.1% women) 
participated, bringing the overall response rate to 78%. The baseline assessments 
included the measurement of anthropometric characteristics, femoral and lumbar spine 
BMD. Subjects were excluded according to the following criteria: age 80 years and over, 
use of a walking aid, use of estrogen or hormone replacement therapy, diuretic, thyroid 

15 hormone or cytostatics, or known diabetes mellitus. After genotyping women with the 
rare VDR haplotypes 4 and 5 (n = 16) were excluded. Anthropometric data, DNA 
samples and genotype data for both loci were finally available in a sample of 
1062 women. Data on incident vertebral fractures were available for a subgroup of 
634 women. The follow-up time was on average 6.5 (SD 0.4) years. 

20 MEASUREMENTS 

At baseline height and weight were measured. BMD (in g/cm 2 ) was measured at 
the femoral neck and lumbar spine by dual energy X-ray absorptiometry (Lunar DPX-L 
densitometer), as reported earlier (2). Body mass index (BMI) was computed as weight 
in kilograms divided by height in square meters (kg/m 2 ). Age at menopause was assessed 

25 by questionnaire. Dietary intakes for calcium (mg/day) and vitamin D (mg/day) were 
assessed by food frequency questionnaire and adjusted for energy intake. Both at 
baseline, between 1990 and 1993 and at the follow-up visit, between 1997-1999, 
radiographs of the spine from the fourth thoracic to the fifth lumbar vertebrae were taken. 
All follow-up photos were analyzed for the presence of vertebral fractures by the 

30 McCloskey/Kanis method (3). The occurrence of non-vertebral fractures was recorded, 
confirmed and classified by a physician. In our study population 85 incident vertebral 
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fractures were captured. 131 subjects had a non- vertebral fracture of which 6 had a 
prevalent vertebral fracture. 

DETERMINATION OF VDR AND ERa GENOTYPES 

For genotyping we determined haplotypes of the BsmI, Apal, and Taql restriction 
5 fragment length polymorphisms (RFLPs) at the 3' end of the VDR gene and haplotypes of 
the PvwII- and Xbal RFLPs in the first intron of the ERa gene by direct molecular 
haplotyping methods as described previously (4). Three frequent VDR haplotypes were 
discerned, encoded 1 (baT), 2 (BAt), and 3 (bAT) (Figure 1). The less frequent 
haplotypes 4 and 5 were excluded from the analysis (n=16). Women carrying these 

10 genotypes represent 1.5% of the population. For direct molecular haplotyping of the 
Pvull and Xbal RFLPs a 346 bp PCR fragment was generated by a forward primer (SEQ 
ID NO:3: 5'-GATATCCAGGGT TATGTGGC A-3 ') and a reverse primer (SEQ ID NO:4: 
S'-AGGTGTTGCCTATTATATTAACCTTGA-S') in a reaction mixture of 10 ^L 
containing 20 ng of genomic DNA, 50 mM KC1, 10 mM Tris-HCl (pH 8.3), 1.5 mM 

15 MgCb, 0.2 mM deoxy-NTP, 2 pM of each primer, and 0.2 U Super Taq polymerase (HT 
Biotechnology Ltd., Cambridge, UK). The reactions were performed in 384-well format 
in a thermocycler (MJ-tetrad) with a cycling protocol of 94, 60 and 72 C for 45 seconds 
each for 30 cycles. Ten microliters of PCR product were digested by simultaneous 
addition of 5 \\\ of digestion mixture containing 5 U Pvull, 7 U Xbal restriction enzyme 

20 (MBI Fermentas) and L5 \xL of ReactBuffer 2 (Life Technologies, Breda, The 
Netherlands) and incubating for 90 minutes at 37°C. The digestion products were 
analyzed by electrophoresis in a 3% agarose gel in 0.5xTBE (lxTBE = 89 mM Tris, 
89 mM boric acid, 2 mM Na 2 EDTA) for 80 minutes at 125 Volts. Separation patterns 
were documented with a digital camera (DC 120, Kodak Company, Rochester, NY, 

25 U.S.A.) under UV illumination (302 nm). Three ERa haplotype alleles were identified, 
encoded 1 (px/T-A), 2 (PX/C-G) and 3 (Px/C-A) combining to six genotypes 11, 12, 13, 
22, 23 and 33 (Figure 2). We did not observe the fourth possible haplotype (pX; - 
397intlT and -351intlG) in our population. 

STATISTICAL ANALYSIS 
30 Differences in mean age at baseline between the study group and the Rotterdam 

Study were evaluated by means of analysis of variance (ANOVA). All other differences 
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in baseline characteristics were compared by analysis of covariance (ANCOVA) testing 
with age to adjust for possible confounding effects. Differences in baseline 
characteristics between the different genotype groups of the ERa gene were compared as 
follows. We grouped subjects by allele copy number (0,1,2) for the haplotype alleles of 
5 interest. We allowed for three possible models to explain differences between groups, 
i.e., an allele dose effect, a dominant effect or a recessive effect. Allele dose was defined 
as the number of copies of a certain allele in the genotype. In case of a consistent trend 
reflected as an allele dose effect, we performed a (multiple) linear or logistic regression 
analysis to quantify the association. In case of a dominant or recessive effect of the test- 

10 allele, ANOVA and ANCOVA tests were performed. For dominant effects we compared 
test-allele carriers versus non-carriers while for recessive effects, subjects homozygous 
for the test allele were compared to heterozygous carriers and non-carriers. 

Odds ratios (ORs) with 95% confidence intervals (95% CI) were calculated by 
(multiple) logistic regression analyses to estimate the relative risk of fractures at baseline 

15 by genotypes of the risk allele, with no copies of the risk allele as the reference group. 
First, we calculated crude odds ratios, and, secondly, we adjusted for potentially 
confounding factors (age, BMI, BMD, and age at menopause). We Used SPSS 
version 9.0 (SPSS Inc., Chicago, USA) for all our analyses. 
RESULTS 

20 BASELINE CHARACTERISTICS 

The study population (n=1062) was on average 67.0 (SD 6.9) years old, had a 
BMI of 26.1 (SD 3.7) kg/m 2 , and age at menopause at 48.7 (SD 4.9) years. Dietary 
calcium and vitamin D intake were on average 1093 (SD 326) mg/day and 1.96 
(SD 1.15) mg/day, respectively. Lumbar spine BMD was on average 1.01 (SD0.17) 
25 g/cm 2 and femoral neck BMD was on average 0.81 (SD 0.12) g/cm 2 , respectively. All 
women were independently living. 

Table 1 shows allele and genotype frequencies for ERa and VDR polymorphisms. 
The genotype distribution was found to be in Hardy Weinberg equilibrium. 
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TABLE 1 Genotype and allele frequencies of ERa and VDR polymorphisms in the 
study population. 



ERa genotype 


Number (%) 


VDR genotype 


Number (%) 


11 


297 (28.0) 


11 


271 (25.5) 


12 


409 (38.5) 


12 


401 (37.8) 


13 


124 (11.7) 


13 


1 05 (9 9*) 


zz 


no ni a\ 
13o (13. U) 


ZZ 


loi (17.2) 


23 


8? f7 7^ 
oz \i *i ) 




8Q (R 4\ 


33 


12 n n 

iz v^.i; 




1 3 n 2^ 


Total 


1062 (100) 


Total 


1062 (100) 


p vdiuc n w i_/ 


U.o 1 


^/ ValUC XiWJj/ 


n 1 1 


ERa haplotype 


Number (%) 


VDR haplotype 


Number (%) 


1 


1127 (53.1) 


1 


1048 (49.3) 


2 


767 (36.1) 


2 


856 (40.3) 


3 


230 (10.8) 


3 


220 (10.4) 


Total 


2124(100) 


Total 


2124(100) 



When we analyzed for known risk factors for osteoporosis by ERa and VDR 
genotypes, no differences were shown apart from ERa haplotype 1 (px) that appeared to 
be dose-dependently associated with later onset of menopause, as we have published 
earlier (Table 2). Similar data were found for the subgroup of 634 women participating 
in the analysis for vertebral fractures (data not shown). 
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TABLE 2 Characteristics of 1062 postmenopausal women according to ERa 
haplotype 1 (A) and VDR haplotype 1 (B). 

A. 



ERa haplotype l b 





Reference 


Heterozygotes 


Homozygotes 


Characteristic" 


(n = 232) 


(n = 533) 


(n = 297) 


Age (year) 


67.6 ±7.1 


66.9 ± 6.9 


66.8 ± 6.9 


BMI (kg/m 2 ) 


26.2 ±3.8 


26.1 ± 3.5 


26.1 ±4.0 


Age at menopause (year) 


47.9 ±5.1 


48.7 ±5.0 


49.2 ± 4.6 C 


Dietary calcium intake (mg/day) 


1098 ±364 


1097 ±319 


1081 ±306 


Dietary vitamin D intake (mg/day) 
B. 


2.04 ± 1.32 


1.96 ± 1.08 


1.89 ± 1.12 


VDR haplotype l b 




Reference 


Heterozygotes 


Homozygotes 


Characteristic" 


(n = 285) 


(n = 506) 


(n = 277; 


Age (year) 


66.6 ± 6.6 


67.3 ± 7.0 


67.0 ±6.8 


BMI (kg/m 2 ) 


26.0 ±3.4 


26.2 ± 3.7 


26.3 ±4.1 


Age at menopause (year) 


48.9 ±4.8 


48.3 ±5.1 


49.0 ± 4.7 


Dietary calcium intake (mg/day) 


1078 ±334 


1104 ±331 


1086 ±306 


Dietary vitamin D intake (mg/day) 


2.03 ± 1.26 


1.91 ± 1.08 


1.97± 1.17 



a Data shown are means ± SD. 

b Reference includes ERa or VDR genotypes 22, 23, and 33; heterozygotes include 12 
and 13; 

homozygotes include 1 1 
10 0 p = 0.02 (allele-dose association, tested by linear regression analysis) 

ASSOCIATION OF ERa AND VDR WITH BMD 

In Table 3 women are grouped according to carrier status for the ERa and VDR 
haplotypes as homozygous carriers (consisting of genotype 11) and heterozygous carriers 
(including the genotypes 12 and 13) of the ERa haplotype 1 (px) and VDR haplotype 1 
15 (baT), respectively, and women not carrying these haplotypes (reference group, including 
genotypes 22, 23 and 33). 
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ERa haplotype 1 was dose-dependently associated with decreased lumbar spine 
BMD corresponding with 0. 1 SD per copy ERa haplotype 1 (Table 3 A, first data 
column). No association was found with femoral neck BMD (Table 3B, first data 
column). ERa haplotype 2 (PX) was associated with increased lumbar BMD, 
5 corresponding with 0.1 SD per copy ERa haplotype 2 (data not shown). These 
associations did not change after adjustment for potential confounders such as age, BMI, 
and age at menopause. No associations were found between ERa haplotype 3 (Px) and 
lumbar spine or femoral neck BMD (data not shown). Based on these data ERa 
haplotype 1 (px) was considered as risk allele. In the sample of 634 women, in whom 
10 data on incident vertebral fractures were available, the association between ERa 
haplotype 1 (px) and lumbar spine BMD showed a similar trend (p = 0.1 1). 

Based on our previous analyses (5) we selected VDR haplotype 1 (baT) as risk 
allele. In the present study no association between VDR haplotype 1 (baT) and lumbar 
spine or femoral neck BMD was observed (Table 3 A and B, first data rows). 

1 5 INTERACTION OF ERa AND VDR GENOTYPES WITH RESPECT TO BMD 

When the association of ERa haplotype 1 (px) with BMD was analyzed according 
to the carrier status for VDR haplotype 1 (baT), there was a significant allele-dose effect 
of ERa haplotype 1 (px) being associated with decreased lumbar spine BMD only for 
women homozygous for VDR haplotype 1 (baT) (Table 3A, fourth data column; p < 

20 0.001). This association was not influenced by age, BMI, and age at menopause. When 
age, BMI, ERa genotype and VDR genotype were taken together in a multivariate 
regression model, there appeared to be a borderline significant interaction between ERa 
haplotype 1 (px) and VDR haplotype 1 (baT) (p = 0.09 for the interaction term). In the 
subgroup of 634 women, in which data on incident vertebral fractures were available, 

25 similar associations were found (data not shown). No interaction between ERa and VDR 
genotypes was found for femoral neck BMD (Table 3B; p - 0.13 for the interaction 
term). 

ASSOCIATION OF ERa AND VDR WITH FRACTURE 
Analysis of the distribution of fractures in women according to the ERa 
30 genotyperevealed an overrepresentation of vertebral fractures in women carrying the ERa 
haplotype 1 (Table 4). 
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TABLE 4 Number of women with vertebral fractures according to ERa genotype. 



ERa genotype 


No. with fracture/total no.(%) 


11 


40/187 (21.4) 


12 


26/236(11.0) 


13 


10/71 (14.1) 


22 


3/79 (3.8) 


23 


6/53 (11.3) 




0/8 (0) 




19.2 


/? value 


0.0002 



Figure 2 shows separately the distribution of vertebral fractures according to the 
ERa haplotype 1 (px) and VDR haplotype 1 (baT) status. Vertebral fractures were 
overrepresented in women carrying ERa haplotype 1 (px). This association appeared to 
5 be dose dependent with 6.4% in non-carriers of ERa haplotype 1 (px), 12% vertebral 
fractures (OR 1.9, 95% CI 0.9-4.1) in women heterozygous for ERa haplotype 1 and 21 
% vertebral fractures (OR 3.9, 95% CI 1.7-8.2) in women homozygous for ERa 
haplotype 1 (px). For women carrying ERa haplotype 2 (PX) there was an allele dose 
association with decreased vertebral fracture risk (p < 0.001) while for ERa haplotype 3 
10 (Px) no differences were observed (p = 0.53) (data not shown). 

In a previous study VDR haplotype 1 (baT) was found to be associated with 
increased fracture risk (5). When women were grouped by VDR haplotype 1 (baT) 
genotype also an allele dose association was observed (Figure 3). Non-carriers of VDR 
haplotype 1 had 11% vertebral fractures, women heterozygous for VDR haplotype 1 had 
15 13% vertebral fractures (OR 1.3, 95% CI 0.7-2.3), while women homozygous for VDR 
haplotype 1 had 18% fractures (OR 1.9, 95% CI 1.0-3.7). VDR haplotypes 2 (BAt) and 3 
(bAT) were not associated with vertebral fracture risk (data not shown). 

No differences between the various ERa and VDR haplotypes were found when 
the risk for incident non-vertebral fractures was analyzed (data not shown). 

20 INTERACTION OF ERa AND VDR WITH RESPECT TO FRACTURE RISK 

When we further stratified by VDR haplotype 1 (baT) genotype we observed the 
ERa haplotype 1 (px) association to be modified by VDR haplotype 1 (baT) genotype 
(Figure 4). Significant ERa haplotype 1 (px) genotype dependent differences were only 
observed in women homozygous for VDR haplotype 1 (baT). Logistic regression 
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analysis showed that, compared to the 'double reference* group, within the group of VDR 
genotype [1,1] women have a 2-fold (95% CI 0.5-7.9) and 10-fold (95% CI 2.7-38) 
increased risk for vertebral fractures when being heterozygous or homozygous for ERa 
haplotype 1 (px), respectively. In non-carriers and heterozygous carriers of VDR 
5 haplotype 1 (baT) no significant ERa haplotype 1 (px) genotype dependent differences 
were observed. When age, BMI, ERa genotype, and VDR genotype were taken together 
in a multivariate regression model, there appeared to be a significant interaction between 
ERa haplotype 1 (px) and VDR haplotype 1 (baT) (p = 0.01 for the interaction term). 
After adjustment for lumbar spine BMD and age at menopause the results did not change 
1 0 (data not shown) . 

DISCUSSION 

The data presented herein demonstrates for the first time interaction of 
polymorphisms in the VDR and ERa gene in relation to the risk of incident vertebral 
fracture risk. Women homozygous for both the VDR haplotype 1 (baT) and ERa 

15 haplotype 1 (px) had a 10-times higher vertebral fracture risk than non-carriers and a 3-4 
times higher risk than carriers of either one of the risk haplotypes. 

So far most association studies focused on single genes. Two meta-analyses 
showed a weak association of VDR genotypes with BMD, which supported our own 
findings in a sample of 2000 men and women from the Rotterdam study. A recent meta- 

20 analysis showed an association between ERa genotypes and lumbar and femoral BMD 
(6). Most genetic association studies for osteoporosis have been performed with BMD as 
endpoint while the clinically more relevant endpoint of osteoporosis is fracture. A 
limited number of studies have yet been able to address the association of specific gene 
polymorphisms with fractures. Previously we have shown that VDR haplotype 1 (baT) is 

25 the risk allele for osteoarthritis, vertebral and nonvertebral fractures (5, 7). Ioannidis 
et al. showed in a meta-analysis that the Xba\ polymorphism and not the Pvull 
polymorphism was associated with increased combined risk for vertebral and non- 
vertebral fractures (6). In the present study direct haplotyping methods were used to 
increase genetic resolution. An association of ERa haplotype 1 (px) with lumbar spine 

30 BMD and vertebral fracture risk was demonstrated. Also at the lumbar spine a 
synergistic interaction between ERa and VDR genotype for BMD and fractures was 
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detected. No interaction effect between ERa and VDR genotypes was found for femoral 
neck BMD and non- vertebral fracture risk. This is in line with previous data, which show 
a higher response to hormonal replacement therapy at the lumbar spine in contrast to the 
femoral neck. The ERa effect may be more pronounced in the spine, which contains 
5 more trabecular bone resulting in a higher rate of bone turnover compared to cortical 
bone, as present for example in the femoral neck. 

We and others previously observed that ERa genotype is associated with 
differences in age at menarche and age at menopause. However, in our current analyses 
age of menopause did not influence the interaction we observed. This suggests that 

10 differences in the age of menopause are small and do not explain the interaction. 
However, because of the relatively small effect, such influences might only be observed 
in studies of sufficient power. 

An interesting observation was that the association of both ERa and VDR 
genotypes with vertebral fracture incidence was independent of BMD. This indicates the 

15 significance of other bone characteristics for the risk of fracture. But it also pointed to 
the involvement of ERa and VDR genes in pathways (e.g., bone matrix synthesis and 
bone turnover) other than those directly reflected in BMD, and which also determine 
strength of bone and thereby fracture risk. For example, estrogen deficiency may 
increase the numbers of remodelling sites and deeper resorption lacunae and therefore 

20 produce loss of connectivity and increased risk for fractures. 

A limitation of the present study may be health selection bias. However, 
genotype and allele frequencies are similar to those observed in other Caucasian study 
populations and so "Health" (apart from the risk for fractures) seems not to be genotype 
dependent and, therefore, we do not expect this to influence the results. Furthermore, 

25 potential selection bias was avoided by deriving cases and non-cases from the same 
source population. 

An aspect that should be realised is that the polymorphisms in the ERa and VDR 
are anonymous. There is no direct known functional consequence for the ERa and VDR 
protein. Therefore, when association is found it is assumed that allele(s) of these single 
30 nucleotide polymorphisms are in linkage disequilibrium with one or more of the truly 
functional polymorphisms elsewhere in the gene. These functional polymorphisms could 
alter VDR or ERa protein structure or might affect the activity of the VDR and ERa 5' 
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promoter and 3 ? UTR, leading to the expression of altered quantities of VDR47 and ERa 
proteins under physiologic conditions. Differential transcriptional activity of the VDR 
and ERa receptor proteins could then preferentially modulate subsets of target genes in 
vitamin D and estrogen responsive pathways. 
5 Although the mechanism(s) for the gene-gene interaction we observe is so far 

unknown, it is conceivable from a physiological point of view. l,25-(OH)2D3 is an 
important factor in estrogen biosynthesis and might thus influence local equilibrium 
between estrogens and androgens. Furthermore l,25-(OH)2D3 regulates ER expression in 
osteoblast-like cells. In this way l,25-(OH)2D3 might regulate the effect of E2 on bone 

10 metabolism. In vitro and in vivo studies have shown that several biological responses to 
treatment with vitamin D, such as intestinal calcium absorption and osteocalcin 
production, are VDR genotype dependent. If l,25-(OH) 2 D 3 influences the effect of E2 on 
bone metabolism, this effect might also be VDR genotype dependent. 

On the other hand, E2 influences vitamin D metabolism and VDR expression. Sex 

15 hormone replacement therapy increases total and free serum l,25-(OH) 2 D3 levels. In 
human fetal osteoblasts E2 upregulates VDR expression. Also in rat duodenal mucosa E 2 
increases VDR expression and bioresponse. In this way E2 might influence vitamin D 
regulated processes, like intestinal calcium absorption and osteocalcin production in 
bone. Several studies have demonstrated that the response to hormonal replacement 

20 therapy is ERa genotype dependent. Therefore, the effect of estrogen replacement on 
vitamin D regulated processes might also be ERa genotype dependent. 

In conclusion, the present study shows an interlocus interaction in relation to 
BMD and fractures between two important candidate genes in osteoporosis. Recently, we 
also demonstrated an interaction between VDR and another candidate gene, the COLIA1 

25 gene, with respect to fracture risk (5). Together, these findings underscore the multigenic 
character of osteoporosis and the importance of the contribution of gene interactions in 
determining fracture risk. At the same time our findings highlight the necessity of large 
(multicenter) studies to achieve sufficient statistical power to further elucidate the 
complex, multigenic character of osteoporosis. 

30 While the preferred embodiment of the invention has been illustrated and 

described, it will be appreciated that various changes can be made therein without 
departing from the spirit and scope of the invention. 
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